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Abstract 



The preheating after the inflation, which can be interpreted as particle creation in the 
oscillating inflaton background, provides a state far from thermal equilibrium. We extend 
the field theoretical treatment of the preheating by Linde et al. to the case of multi- 
component complex scalars to show that some charges are created in this process, if C 

* i 

and CP are violated. A new possibility of baryogenesis based on this mechanism is also 
H \ discussed. 
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1 Introduction 



In the scenario of inflationary cosmology, which can explain the unanswered problems left 
by the standard big bang scenario, all the particles at present are thought to be created in 
the reheating process following the exponentially inflating de Sitter phase. That process 
is viewed as the decay of the inflaton field, which begins to oscillate around the minimum 
of the inflaton potential. The inflaton itself is treated as a coherently oscillating classical 
field, and its oscillation is damped by its decay, whose decay products fill the universe to 
realize the hot universe. A several yeas ago, it was pointed out by Linde et al.|I[ that 
light bosons could be produced more drastically than the perturbative reheating, when 
some conditions are satisfied on their coupling to the inflaton and the amplitudes and 
frequency of the inflaton oscillation. This process is called preheating, to distinguish from 
the perturbative reheating. 

The preheating is nothing but particle creation in the oscillating background. Al- 
though particle creation or annihilation generally occurs in an time-dependent back- 
ground, the number of created particles grows exponentially when the background is 
periodic in time. This is a kind of parametric resonance. This process of particle creation 
is far from equilibrium, and the state after that has large quantum fluctuation composed 
of low frequency modes, whose relaxation process to the equilibrium is also out of equi- 
librium. Hence there are chances for the matter-antimatter asymmetry of the universe 
to be generated in these eras, if C and CP are violated during these processes. In fact, 
the preheating has been argued in connection to generation of the baryon asymmetry of 
the universe (BAU)0, focusing on the fact that the low-frequency quantum fluctuations 
of the created particles are so large that the sphaleron processes may be enhanced. Be- 
sides this possibility, it is natural to think that the baryon number or other charges were 
generated during the preheating when all the particles in the universe were created. We 
investigate this possibility to generate a charge asymmetry in the presence of C and CP 
violation in the oscillating background. 

We study the particle creation when C and CP is violated in the interaction of quan- 
tum scalar fields with the background. We extend the analysis by Linde et al.|| to a 
system of multicomponent complex scalar fields. The interaction with the oscillating 
background is approximated by the successive scattering with some potential, while the 
effects of CP violation are evaluated perturbatively. In ||, the number density of the 
created particles during the preheating was expressed in terms of the wave function in the 
inflaton background, and the results were reproduced by the picture of successive scat- 
tering approximation. It is, however, difficult to distinguish the positive-frequency wave 
from the negative-frequency wave in the oscillating background, so that one cannot define 
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a charge, to which the contributions from the antiparticles have opposite sign relative to 
the particles, in terms of wave functions. Then the treatment based on the successive 
Bogoliubov transformations will be useful. In Section 2, we generalize the field theory of 
a real scalar in the oscillating background studied in || to the case of multicomponent 
complex scalars and formulate the descent equations for the Bogoliubov coefficients. We 
apply the method to a simple toy model and find some charges to be generated, in Section 
3. The final section is devoted to discussions. 

2 Field Theoretical Treatment 

Consider a system of n complex scalar fields interacting with the inflaton through the 
effective lagrangian 

£ = d^Xa - m 2 a (t)xlXa - xlV ah (t) Xb - \ (XaW ab (t) Xb + c.c.) , (2.1) 

where m 2 a {t) = g 2 (p 2 {t) and a,b = 1, 2, • • • , n. The coefficients V a b(t) and W a b(t) represent 
the collective effects of the oscillating background and interactions with other fields. We 
treat their effects as perturbation. For simplicity, we consider the case of the undamped 
background oscillation, with the parameters in the broad resonance regime. That is, we 
assume that the inflaton field is given by 

= $ sinmt, (2.2) 

where m is the inflaton mass and $ is the constant amplitude, and neglect the redshift 
of frequency of each mode. In this case, each mode of the unperturbed wave functions 
satisfies the mode equation 

(d 2 + k 2 + g 2 y sin 2 mt) X ak(t) = 0. (2.3) 

It is known that wave functions in any periodic potential belong to either stability or 
instability bands[|]]. The solutions to ( |2.3| ) are given by the Mathieu functions, which 
have this property. Among the solutions, those in the instability bands increase their 
particle numbers exponentially. In particular, for q a = ^-g- 3> 1, some modes are in the 
broad instability bands so that they continue to grow for a long time, even if we take the 
effects of the expanding universe into account. Further, the frequencies of the modes in 
the broad resonance band are much larger than that of the inflaton oscillation. Then, 
the particle number changes only when the inflaton field vanishes, while the x-fields 
evolve adiabatically for time intervals with \4>{t)\ ^ 0. This observation motivates the 
study of preheating by approximating the oscillating background as a series of successive 
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scattering by some simple potential. In the inflaton field near its zero was replaced 
by the quadratic potential, and the descent equation for the particle number density 
in some mode was derived. The particle number density is given in terms of the wave 
function, which is a solution to ( |2.3|) with any boundary condition. As for a charge 
density, we must distinguish the particle from the antiparticle. This amounts to impose 
the boundary condition defining either positive- or negative-frequency mode. Here we 
express the densities of particle number and the charges in terms of the Bogoliubov 
coefficients which relates the operators before the inflaton begins to oscillate and those in 
the adiabatic time intervals. 



2.1 Descent equations for the Bogoliubov coefficients 

We assume that the particles are generated only in the short time interval near tj = nj/m a , 
where j is a non-negative integer, and that the violation of the charge conservation and 
CP symmetry is effective only within the short interval when m a (t) ~ 0. Hence, for t 
with tj-i <C t <C tj, Xak{t) satisfies the unperturbed equation (|2.3|) . Let f 3 ak {t) be the 
properly normalized positive-frequency solution to (|2.3|) in the adiabatic interval between 
tj-i and tj. Then the \- field operator is expanded as 

Xa(x) = J d 3 k (ai k f ak (t)e ik * + Klfil(t)e- ikx ) . (2.4) 

Without the perturbation, upon scattered by the potential g^(j) 2 (t) at tj, the positive- 



frequency wave f 3 ak (t) is transmitted to a linear combination of f 3 ak (t) and f 3 ak *(t), 
whose coefficients are determined by solving the scattering problem by the potential at 
tj. In the presence of the perturbation V a b and W a b, the scattered wave is composed of 
various components of f 3 ak l {t) and f 3 ak l *{t). Suppose that the positive- frequency wave 
fa k (t) at t < to transfers to a set of the positive- and negative- frequency waves in the j-th 
interval; 

/ o \ ( <j( k (t) + pijf k (t) \ 



f°ak(t) 



< b fL(t) + /&/£(*) 



(2.5) 



while the negative-frequency wave fa k (t) transfers as 

/ o \ / puis) + aiji*(t) \ 



v o / 



(2.6) 
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where a 3 ab ^ a 3 ab and (3 3 ab ^ f3 3 ab because of the CP violation. Then the Bogoliubov 
transformation between the creation and annihilation operators at t < to and those at 
t 3 -_i <t <tj is 



u ak 



a bk a ba ^ u bkl J bai 
a bkPba + ^bk^ba- 



Requiring 



j j] 

a aki a bk' 



KkXi]=^(k-k'), 



fl afc5 Kk' 



(2.7) 



(2.8) 



we find the conditions on the Bogoliubov coefficients: 

«ca«cfe - PLPlb = $ab, 

alaPlb - PL°& = °- 



(2.9) 



In terms of the n-hy-n matrix notation, a J = [a 3 ab j , (3 J 
expressed as 



lj ab ), etc., these conditions are 



cu 3 V - (3 j ^(3 j = a 3] a 3 - (3 j ^(3 3 = 1, /?V - a j] L = 0. 



We can also show that 

oc j a j] - {3 3 f3 j ^ = a 3 a 3] - (3 3 ^ = 1, 



a 



j {3 3 ^ - fta j] = 0. 



(2.10) 



(2.11) 



Now we define the number density of the particles with momentum k in the j-th 
interval as 



i n 

^(o°lE(^ fc + ^ fc )|o ), 

1 a=l 



[2.12) 



where |0°) is defined by a° fc |0°) = b° ak \0°) = and V = 5 3 (k = 0). By use of we 
have 

4 = E {&L + PiaPil) = Tr + fPffi) . (2.13) 



a.b 



If each component of (f> a {x) has the U(l) charge Q a , the density of the generated charge 
in mode k in the j-th interval is given by 



V a=l 



which is written in terms of the Bogoliubov coefficients as 



fk 



Qa [PbaPba PbaPba 



Tr 



Q((3rt[3 j - fit & 



(2.14) 



(2.15) 



where Q in the last line is n-by-n matrix, Q = diag(Qi, Q2, • • ■ , Qn)- 

In the case of n — 1, it is easy to find \ft\ 2 = ft by use of (|2.10|) , so that j{ = 0. 
Hence, we need at least two complex scalars interacting with the oscillating background 
in order to generate charge asymmetry. This fact can be understood as follows: If we 
regard the particles and the charges as the decay products of the oscillating field, at 
least two decay channels are needed to produce nonzero charge since the effect of CP 
violation appears as the interference of the two channels. This is the same situation as 
the GUT-baryogenesis by the heavy boson decay. 

The Bogoliubov coefficients in the (j + l)-th interval are related to those in the j-th 
interval by solving the scattering problem by the potential at tj, just as done in ||. Within 
each adiabatic interval, the positive frequency wave function f ] ak {t) is well approximated 
by 

h ak (t) = 1 e -ifi<#M«) } (2.16) 



/2w a (t) 
where 

u a {t) = ^/k 2 + m 2 a (t) = y/ 'k 2 + g 2 ® 2 sin 2 mt> (2.17) 

up to some constant phase. On the other hand, for t near tj, for which the adiabatic 
approximation is no longer valid, we solve the scattering problem employing some ap- 
proximate form of the potential at tj. We adopt another type of potential which allows to 
define an asymptotic state, that is, uj a (t) ~ uj® = const, for \t — tj\ ^> 1/m, in contrast to 
the quadratic potential rn? a (t) ~ m 2 (t—tj) 2 used in [g]. This choice of the potential is more 
appropriate for the perturbative calculation of the CP- violating scattering presented in 
the next subsection. The general formulas presented in this subsection holds independent 
of the detailed form of the potential. If we use this potential, the adiabatic wave (|2.16|) 
behaves, far from the potential, as 

h ak (t) -> -^ e -^(*-%)-<*-<*, (as t - tj — ±00) (2.18) 



where 6 J ak = Jq dt' u a (t') is the phase accumulated by the moment tj and $+ fc ($^ fc ) is 
the phase difference between f t 3 dt' u a (t') and u®(t — tj) for t ^> tj (t tj). Then the 
positive- frequency wave at t < to, fak(t)i can be approximated by h ak (t), while f ak {t) in 
the j-th adiabatic interval will be expressed as a linear combination of hb k (t) and hl k (t). 

Suppose that we find the solutions to the wave equations (usually, by perturbative 
method) satisfying the boundary condition that only the positive-frequency wave in the 
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a-th component exists at t ^ tf 
( \ / \ 



J_ „-iwg(t-tj) 



V 



t<t, 



Xa*(*) 
/ V Xnk(t) J 



t>ti 



2w° 



•^an /- — 77 ° '-'an rz — 77 ° 



(2.19) 

and similarly the solution satisfying the boundary condition that only the negative- 
frequency wave exists at t <C tj\ 



( 

























\ Xn*(*) / 





/ b i — i p-H'C*-^) + i , L= f>«"°(*-*i) \ 



V 



(2.20) 

Here the coefficients A, B, A and B, which are also kinds of Bogoliubov coefficients, satisfy 
the same relations as a, /3, a and f3: 



A ] A - B ] B = A ] A - B ] B = 1, B ] A - A ] B = 0, 
A4 f - BB^ = AA ] - B& = 1, A& - BA 1 = 0. 

Inverting ( j2.5|) and (|2.6| ) and replacing /° fc (t) with h a k(t), we have 



(2.21) 



and 



( ^ 










fL(t) 






)>(*)-( 




V o ) 




I («*) 


an V 




( ^ 








/?•") , M*) \ 

/ al 










' ao 


v o y 








K k {t) j 

J an ' 



(2.23) 



Since how h a k(t) and h* ak {t) appearing in each element in the right-hand sides of ( |2.22j ) and 
Q2.231 ) are scattered by the potential at tj are determined by ( 2.19Q and (|2.20|) , respectively, 



'2.22) 
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the 6-th component of the right-hand side of ( j2.22| ) becomes, after the scattering (t ^> tj), 
(art) e- l ^ + ^A cb -(^\ e^+^BJ e^+^ h bk (t) 

V / ac V / ac J 

W e l ^ + ^A cb -(art) e-^+^BJe-^+^hUt) (2.24) 

: V / ac -I 



which should be identified with (V +lt )^ h bk (t) - {f3 j+1] ) nh Kbit), and the b~ th component 



of the right-hand-side of ( |2.23|) becomes 
'(art) 



e^ + ^B cb - (3rt) e-^+^A 



■cb 



((3rt) e -^ ck +^ k ) Bcb _ (<pf) e^^Aj e^+O h* bk (t). (2.25) 

V / ac V / ac J 

which should be identified with (a j+1 ^ ^ h bk (t) — h bk (t). If we denote the n-hy-n 

matrix of the phase factors as 

U ± = diag (e'^ + <\ • • • , e-^ + ^) , 

we obtain the descent equations for the Bogoliubov coefficients 



(2.26) 



a 



= u+ (3 j - &U-a? 

a j+1 = Ul(A^U_a j -tfulp* 
(3 j+1 = U{(A ] U-ft -B ] Ula 3 ) . (2.27) 

As a check, one can show that these (j + l)-th Bogoliubov coefficients satisfies the same 
relations as Q2.1CQ and (|2 . 1 l|j . 



2.2 Perturbative calculation of the scattering matrices 

Now we calculate the scattering matrices A, B, A and B by the perturbative method. 
The wave equations to be solved are written as 



(d 2 t +e + ml(t))l 



(V ab (t) W* ab {t) 

\w ab (t) v: b (t) 



( Xbk{t) 

\xt k (t) 



0. 



For convenience, let us introduce the following symbols: 



X k {t) 



Xakit) 
X ak (t) 



A = (% + tf + ml{tj)l= 1, 

v ab {t) W* ab {t) 



A! 



W ab (t) V: b (t) 



(2.28) 



(2.29) 
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Then a perturbative solution to the equation 



(A + Ai)X k (t) — 0, (2.30) 
to the first order with respect to A± is given by 

X k (t)=xl°\t)-A 1 A 1 xl°\t), (2.31) 
where X ( k °\t) is a solution to 

A xi 0) (t) = 0. (2.32) 

In order to evaluate the scattering matrices A and £>, we must prepare Xk(t) such that 
only Xak(t) is the positive-frequency wave and the other components vanish in the remote 
past. We solve the scattering problem with this initial condition for the approximate 
potential which admits the asymptotic states. Let f a k(f) be the wave function satisfying 

fe 2 +^(t)U fc (t)=0, and f ak {t) *=i=S° -^=e~*#. (2.33) 



Suppose that we solve this unperturbed equation and find the expression as 

fauit) = a a fT(t) + (3JT(t), (2.34) 

where satisfies 

(^ 2 + ^))/rW = 0, and /a 7(*)^ 1 e"^. (2.35) 



According to the general theory of distorted- wave Born approximation 0, the Green func- 
tion G a (t,t') satisfying 

(d^ + ul(t))G a (t,t') = 6(t-t'), (2.36) 

with the boundary condition such that there is only positive-frequency wave in the remote 
past is given by 

(~i (f 4-1 \ _ J w fak{t)fak\t )' \t ^ ^ ) /'o q7\ 

~ i (*'<*) { } 

where w is the Wronskian defined by 

W = Umfakit) - flkiWakit), (2.38) 

which is time-independent and w = i in our case. Then our Aq 1 is 

A 1 = diag (G 1 , ■ ■ ■ , G n , Gl, ■ ■ ■ , G*) . (2.39) 
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As the unperturbed wave X^ k \t), we take 



M(o)aa _ ( c fk( f ) 



/Calflk(t)\ 
<if* k (t) 



\ 



(2.40) 



/ 



and chose the set of parameters c ab such that 



X [ k a] (t) = X l ° m (t) - A^AiXp^t) *: 



/ \ 

fak(t) 



(2.41) 



V o J 

In terms of the n-component wave function, the perturbed wave becomes 

xt\t) = cf k (t) - GV{cf k ){t) - GW\c*r k ){t). (2.42) 
Requiring the initial condition Q2.41 ), we find that c a t should be taken to satisfy 

poo poo 

c ab + ic ac / dt' v bc (t')K k (t')f ck (t') + ic* ac / dt' w; c (t')f; k (t')r ck (t') = s ab . (2.43) 



Once is determined, we obtain the scattering matrices A and B by extracting the 
coefficients of the positive- and negative-frequency waves, respectively, from ( [2.42 ) in the 
remote future: 



A ab = c ab a b + i/3; dt'f bk (t')[V bc (t') C(u: f ck {t') + WUt')clj: k (t')}, 

J — oo 

POO 

B ab = c ab (3 b + ia* b ^/^[^(Ocac/dkCO+W'teCOc^CO]- ( 2 - 44 ) 



Similarly, by preparing the negative- frequency wave in the remote past, we can find the 
other set of the scattering matrices A and B. They are given by 



A ab = 5 ab a* b -t(3 b dt'n k {t') [V bc (t , )c ac f ck (t , ) + W^{t')cUck{t')\ ■ 

J — oo 

/oo 
dt'fuf) [v bc {t')cacf: k {t') + w: c {t')cij ck {t')} 
-oo 

where c ab must satisfy 

/•oo roc 
Cab - ICac \ dt' V bc (t')f bk (t')P ck (t') - l~C* ac / dt' Wfc(0/wfc(0/«*(0 = fo- 



(2.45) 



(2.46) 
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In the case of W a b(t) = 0, for which the global U(l) symmetry exists corresponding to 
the rotation of all Xa{x) with a common phase, the scattering matrices are simply given 
by 

A = N- 1 (a + Kf3 ] ) , B = N^ 1 (j3 + Ka^ , 

A = N- 1 ^ + K/3), B = N- 1 (/3 f + Ka) , (2.47) 

where a and (3 are diagonal matrices defined by 

a = diag(ai,a 2 ,---,a n ), P = diag (ft, fa, ■ ■ ■ , P n ) , (2.48) 
and the matrices N, N, K and K are given by 

/oo 
dtV ba {t)f ak (t)fZ k {t), 
-co 

/oo 
dtV ba (t)f: k (t)f bk (t), (2.49) 
-oo 

/oo 
dtV ba (t)f ak (t)f bk (t), 
-oo 

/oo 
dtV ba {t)r ak {t)r bk {t). (2.50) 
-oo 

Because of the hermiticity of the matrix V ab (t), = K and T ab = Jf^ c?t V ba {t)f ak {t)f b \{t) ■ 
T b * a , the matrix can be regarded as the unitary matrix exp(iT). 

3 Example 

Once the model of the inflaton interacting with at least two complex scalars is specified, 
one can evaluate the 'effective potential' V ab (t) and W ab (t), which may be induced directly 
by coupling with the inflaton to the complex scalars, or through some loop effects including 
the scalars and any other fields. Here, without proposing any specific model, we illustrate 
calculation of the generated charges and particle numbers, assuming some simple form 
of the CP- violating potentials. For simplicity, we concentrate on the case of W ab (t) = 
and n = 2 and assume the degeneracy of the scalar masses (mi = 777-2) and the property 
of the potential V\\ = V22, which may result from some discrete symmetry. Under these 
assumptions, 

fik{t) = f2 k {t) = /*(*), ai = a 2 = a, p 1= p 2 = p. (3.1) 
Then the matrices in fl2.49|) and (|2.50| ) are reduced to 
K -- 



%K X iK 2 {l-ie K )\ K = K ] 

iK 2 (l + ie K ) iK x 



N 



( 1 + m ^ iN 2 (l-ie N )\ ^_ t 
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where 

/oo roo 

dtV n (t)\f k (t)\ 2 , N 2 = / dtReV l2 (t)\f k (t)\ 2 , 
-oo J — oo 

e ^ = W r dtlmV 12 (t)\f k (t)\ 2 , (3.3) 

iV2 J—oo 

/OO /"CO 

dtVuWZit), K 2 = / diReV 12 (i) / 2 (t), 
-oo J—oo 



-oo 

roo 



e K =^r I dt\mV 12 {t)f k {t). (3.4) 



K- 



Once the explicit forms of /&(£), Vii(t) and Vi 2 {t) are provided, we can calculate all the 
quantities needed to follow the descent equations ( |2.27| ), so that the time evolution of the 
densities of particle number and charges. 

First of all, we prepare the unperturbed wave which is appropriate for the perturbative 
calculation of the scattering matrices. For this purpose, we substitute the oscillating 
potential near its zeros with the asymptotically flat potential as 

$ 2 sin 2 mt ~ 2$ 2 tanh 2 ^ ^y=~^ j ■ ( 3 - 5 ) 

This is better than the simple quadratic potential adopted in 0, in that the potential 
of ( p.5|) matches to the original oscillating potential up to the forth order of m(t — tj). 
Introducing the dimensionless parameters 



k I m 1 

r = K(t-t 3 ), «=-, l=^— = w „ (3.6) 

where k* = g$m = y/^mq 1 ^ and q = the wave equation is expressed as 



d 2 h{r) 
dr 2 



K 2 + ^tanh 2 (£r) 



fk(r) = 0. (3.7) 



We find that the solution which asymptotically approaches the positive-frequency wave 
in the remote past is given, in terms of x = £r, by 



v ' \ 1 + e- 2 *, 

V 2uJ o 



where 



9 ^ k + ^— = 2q" 2 Jl + ^— 



m 

uq = 2el;k* = \/2me, (3.9) 
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VV V/2 V 64g/ ' 2 



K 2 



c = -W /2 f + 4^72+1, (3.10) 

and -F(a, 6, c; z) is the hypergeometric function. By use of its property, we find, in the 
limit of t — t 3 • — >• oo, 

f k {t) — > _^= e — o(^) + A e «ao(^) (3 n) 

V-M) v^o 

where 

T(-2ie+ l)r(-2ie) 

OLk 



T (I - 2*e + 4^1/^1 - r (1 - 2*e - 4^/2^1 - ^ 

r(-2^e + l)r(2ie) . cosh (W^ 1 ~ efc 

Pfc - 



r(l + 4V/ 2 ^4)r(i-4V/ 2 v / ^4) sinh ( 27re ) 

(3.12) 

Among the phases defined in ( 2.18| ), is given by 

0£= P dtu(t)=j9 k , (3.13) 



with 

6 k = 2V2q 1/i £ d^K 2 + 2gV2 sin 2 ~ ( 3 . 14 ) 
while is evaluated by use of the potential ( |3.5|) as 

= 2e log ' - 4^ log fV4(e + V/2) (3.15) 

With these quantities, one can evaluate the density of the particle number in the absence 
of the perturbation. The result reproduces the same as that in ||, for the case of q 1 ^ 2 = 
l/(4£ 2 ) » 1. 

As for the effective potential, we adopt 

V n (t) = V 22 (t) = 2A 1 $ 2 (tanh 2 £r - 1), 

V 12 (t) = 2A 2 $V e{r) (tanh 2 £r- 1), (3.16) 



which have the same t-dependence as the potential in ( |3.7|) and vanish for large \t\. We 
will specify the CP- violating phase 6(r) later when we execute numerical calculation. If 
we denote 

h = \ k = \ (3.17) 
9 9 
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the quantities composing the elements of the matrices ( p.2| ) are given by 



No 



2kg 1 / 2 



oo ] 

dx ■ 



1 + K 2 /(4gV2i .'-x 



^1 + K 2 /(4g 1 /2) 

1 2/ 2 gV 2 
^2 + K 2 /(4gV2) 

2/xg 1 / 2 
^/l + K>/{4qW) 
2kg 1 / 2 



dx 



cosh x 
cos6*(x) 



F [a,b, c; 



1 



cosh 2 x 
sin 9(x 



F ( a, 6, c; 



dx 



cosh x 



F ( a, 6, c; 



1 + e" 2a; 
1 

1 + e~ 2x 
1 



1 + e 



-2x 



dx . 
oo cosh X 



00 , cos#(x 
ax 



1 _|_ K 2 /(Ag l l 2 ) J -°° cos h x 
I 2/ 2 g 1 / 2 



(> + 



He f 1 

~' M F ( a, 6, c; 
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Now we are ready to evaluate the densities of the particle number and charges. The 
remaining quantities we must specify are the values of q, l\ and I2, and the functional form 
of 0(t). All the quantities with mass dimensions are expressed with appropriate powers 
of the inflaton mass m. For definiteness, we adopt 
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6{x) 



1 + e 2 * ' 



(3.19) 



which changes from 8q to as x goes from — 00 to 00. We take g = 200, for which the 
parametric resonance predicted by the Mathieu function occurs for a wide range of the 
momentum k, that is, in the broad resonance. We numerically evaluate the scattering 
matrices ( 2.47|) for a fixed k, and study the time evolution of 
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(3.21) 



where (|3.20|) and (|3.21|) are the charge densities defined with Qi 


1 

the initial condition a 



1 




and Q2 



in ( 2.15Q , respectively, by numerically solving the descent equations ( 2.27 ) with 



a = 1 and (3° — (3° = 0. Because of the U(l) symmetry 
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mentioned above (|2.47|) , j{ k + f 2 k = in any j-th interval. Although the particle number 
grows exponentially for k in one of the resonance band, it ceases to grow due to the back 
reaction after a finite number of inflaton oscillations, as studied in ||. How many times 
the inflaton oscillates before the particle creation ends depends on q and the coupling 
constant g. For simplicity, we use the densities in the tenth interval when we evaluate the 
total densities of the generated particles and charges, since 10 is a typical number of the 
oscillations before the preheating ends. 

Now let us show the numerical results. At first, neglecting the back reaction, we 
studied the time evolution of the densities for 6*o = 10~ 3 , l\ = 0.01 and l 2 = 0.02. Since 
(3k in ( |3.12| ) behaves as \(3k\ — e -7 ™ 2 , the modes with k satisfying k < — 0.56 are 

expected to belong to the first resonance band. The time evolution of ri]~ and jik for 
several k's are shown in Figs. [I] and 0. As expected, the densities for k < 0.5 grow 
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Figure 1: Time evolutions of the densities n k (solid curve) and jik (dashed curve) for 
k = 0.1 and 0.3. 



exponentially, and 



10 4 nl for this choice of the parameters. 



We estimate the total amount of densities of the generated particles and charges by 
integrating n J k and j[ k at j = 10 over k: 



d 3 k n k = 8V2' 

8V27T m 3 g 3/4 / dnn 2 j 
Jo 

POO 

j 2 = 8V27rmV /4 / dKK 2 j 2k . 

Jo 



7rm 3 g 3//4 / dKK 2 n k , 
Jo 

00 



Ji 



(3.22) 



We show the behavior of the integrands at j = 10 as functions of k for several values of #0 
in Fig. |3|. Except for the case of 9 = ir, the generated particle numbers are not affected 
by the CP violation, while the generated charge is almost proportional to 6q. As seen 
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Figure 2: Time evolutions of the densities n k (solid curve) and j\k (dashed curve) for 
k = 0.5 and 0.6. 



#0 


J dn n 2 nk 


/ dn n 2 j lk 


10" 3 


130.5096 


-1.609334 x 10~ 2 


io- 2 


130.5156 


-1.544579 x IO" 1 


io- 1 


131.1163 


-1.537716 


7T 


990.7411 


-50.84228 



Table 1: The integrals appearing in the total densities of particle number and charges 
( 3-22 ) for several values of 9 . 



from the K-dependence of the integrands of ( 3.22 ) in Fig. y, the second resonance band 
exists near n = 0.8 but its contribution is about five orders smaller than that from the 
first resonance band. We performed the calculation up to n = 10 and found several higher 



resonance bands, whose contributions are negligible. The values of the integrals in ( 3.22 
are listed in Table [|. 

We also calculated the total densities for various values of the coupling constants l\ 
and I2 for 6q = 7T. As shown in Table there is no regularities in how the total amount 
of densities depend on the coupling constants. 



h 


/ dn K 2 rtk 


I dK K 2 j lk 


10" 


-i 


1.071055 x 10 5 


2.456463 x 10 2 


10- 


-3 


6.773442 x 10 4 


-2.316711 x 10 3 


10" 


-2 


9.907411 x 10 2 


-5.084228 x 10 1 


10- 


-1 


1.656527 x 10 3 


-2.069592 x 10 2 



Table 2: The integrals appearing in the total densities of particle number and charges 
( 3.22| ) for several values of l\ and I2 = 1l\ with 9q 
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Figure 3: The densities of the particle number n 2 nk and charge —n 2 jik at j = 10 for 
#o = 10~ 3 (solid curve), 10~ 2 (dashed curve), 10 _1 (dotted curve) and n (dashed-dotted 
curve). 

4 Discussions 

Since all the particles at present are thought to be created after the reheating process in 
the inflationary scenario and the process is far from equilibrium, it is natural to consider 
that the asymmetry between matter and antimatter was also generated in that era. We 
formulated charge generation during the preheating, which is viewed as a particle creation 
in the oscillating background. 

If the scale factor of the FRW universe changes from -R pre heating to R th just after 
thermalization, the charge-to-entropy ratio will be given by 

jl _ 8V2ll m 3 g 3 / 4 J °° dK K 2 j lk ( ^preheating \ 3 , . 

S "45" -9*S-l th V K th / 

where g^s counts the degrees of freedom contributing to the entropy and T th is the tem- 
perature just after thermalization. So far, we have not specified what the charge is, but we 
expect it to be conserved after the preheating process. One may think it as the baryon or 
lepton number, regarding the x-held as the scalar partner of the quarks or leptons. If the 
preheating occurs before the electroweak phase transition, the generated lepton number 
is converted to the baryon number through the sphaleron process. When the charge is the 
weak hypercharge, the equilibrium state after thermalization with the hypercharge excess 
will become baryon-rich state, just as predicted in the theory of electroweak baryogenesis. 
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Anyway, our next subject will be to construct a realistic model, including at least two 
complex scalars, which admit the charge generation based on the mechanism studied here 
and to evaluate the baryon asymmetry produced by that model. 

Although we have considered the coherently oscillating field as the inflaton, the for- 
malism may be applicable to the Affleck-Dine scalar, which is coherently rotating in the 
internal space with a constant frequency if we neglect the damping due to the expansion 
of the universe. The AD scalar, which is the classical squark or a slepton field, is thought 
to decay perturbatively to quarks and leptons to leave B or L. This situation is similar to 
that of the conventional reheating process and one may expect that there is also a process 
similar to the preheating in the AD mechanism. 

Acknowledgments 

The authors thank J. Sato for valuable discussions. This work is supported in part by 
Grant-in-Aid for Scientific Research on Priority Areas (Physics of CP violation, No. 12014212) 
from the Ministry of Education, Science, and Culture of Japan. 

References 

[1] L. Kofman and A. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73 (1994) 3195. 

[2] J. Garcia-Bellido, D. Grigoriev, A. Kusenko and M. Shaposhnikov, Phys. Rev. D60 
(1999) 123504. 

[3] L. Kofman and A. Linde and A. A. Starobinsky, Phys. Rev. D56 (1997) 3258. 
[4] L. D. Landau and E. M. Lifshitz, Mechanics, Section 27. 

[5] See, for example, K. Funakubo, A. Kakuto, S. Otsuki, K. Takenaga and F. Toyoda, 
Phys. Rev. D50 (1994) 1105. 



18 



